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Abstract 

The existing derivations of a heavy quark effective theory (HQET) are analyzed 
beyond the next-to-leading order in 1/m. With one exception they are found to be 
incorrect. The problem is a wrong normalization of the heavy quark field in the 
effective theory. We argue that the correct effective theory should be given by a 
Foldy-Wouthuysen type field transformation to all orders in 1/m. The renormal- 
ization of the resulting Lagrangian to order 1/m^ is performed including also effects 
arising through vacuum polarization. Our results for the anomalous dimensions dis- 
agree with the existing ones. Some applications are considered. 
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1 .Introduction 



Interest in the physics of the mesons and baryons containing one heavy quark (c or 
b) has received a strong stimulus through the invention of the heavy quark effective 
theory (HQET) This allows for a systematic description of the symmetries 
appearing in the infinite quark mass limit and of their breaking by a finite quark 
mass. An important outcome of this theory concerns the possibility of accurately 
extracting the Kobayashi-Maskawa matrix element Vet in a model-independent way 
through a study of the weak semileptonic decay B D*eu. This is made possible 
by Luke's theorem p|, which asserts that some of the 1/mc-order corrections to the 
infinite mass limit predictions of the theory for this decay vanish when the c quark 
is produced at rest in the b quark's reference frame. The next corrections come 
thus in at order and their estimation is clearly important for an assessment 

of the reliability of the Vcb determination. The task of evaluating these corrections 
has been recently undertaken by Falk and Neubert , who made use for this of the 
HQET formalism developed in The aim of this paper is to point out that 

these treatments are incorrect, already at order We show in Section 2 that 

the problem consists in a wrong normalization of the field in the effective theory. 
The correct treatment turns out to be the one given by the authors of . However, 
due to the particular strategy of calculation adopted in 0, their results remain 
valid when the correct treatment is being used. In Sect. 3 the matrix elements of 
current operators obtained with the two HQETs are directly compared and shown 
to differ at order A general presentation of the HQET as given in to any 

order in 1/m is offered in Sect. 4. The resulting effective Lagrangian to (9(l/m^) 
is used in Sect. 5 to sum up the leading logarithms of the heavy quark mass of the 
form l/m^log(m) in the various matrix elements of the theory. We compute also 
the quantum corrections appearing to this order from graphs containing one heavy 
quark inside closed loops, the analog of the Euler-Heisenberg Lagrangian relevant 
to the infrared behaviour of QCD. An additional complication which appears when 
light quarks are included is pointed out: light-quark operators are induced through 
renormalization. In Sect. 6 applications of the formalism to the phenomenology of 
heavy hadrons are given. 



2. Normalization of the field 

Three distinct derivations of a heavy quark effective theory (HQET) from QCD have 
been proposed in the literature 0, ||, ^. Two of them 0, |[, although phrased in 
different formalisms, are actually equivalent. The third one gives different results 
from the first two, starting at order 0{l/m'^). The aims of this section are to explain 
the reasons for this disagreement and to show that the last theory is the correct one. 
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We start by first reviewing the results of the usual derivation of the HQET as 
given in [^]. The relation between the heavy quark field in QCD, Q{x), and that in 
the effective theory, h{x), is up to (9(l/m^) 



Qix) 



—Hp 



(1) 



Here is the velocity of the heavy quark, = d^ + igA'l^t"' and D'^ = D'^ — v^v ■ D. 
The heavy quark field h{x) satisfies yTi = h. The fiavour- conserving vector current 
for quarks with the same velocity can be written by making use of (1) as 



h,h + — 



+ 



1 ^ 



Am? 



K-i^v Dp^ + P^ 7a. Pi. + P±v- D i^,)h . 



The effective theory Lagrangian is 



C = h 



... D P] 



Am? 



P^ V- Dp_ 



h. 



(2) 



(3) 



Let us consider the number operator for heavy quarks at rest {v^ = (1,0,0,0), 
7o/i(x) = h{x)): 



N = / d^x{Q-foQ){x) = / d^xh{x 



1 + 



4m^ 



Pi 



h{x) 



(4) 



Its expectation value in a hadronic state containing the heavy quark Q, also at rest, 
is 



1 



2m 



M 



{Mqcd{v)\ d?x{Q^oQ){x)\MQCD{v)) 



(5) 



2m 



M 



-{Mhqet{v)\ / d?xhH{x 



Air? 



Pi 



hH{x)\MHQET{v)) . 



In this relation the field hnix) is in the Heisenberg picture and its evolution is dic- 
tated by the full HQET Lagrangian (3). \Mqcd{v)) is an eigenstate of the QCD 
Lagrangian and \Mhqet{v)) is the corresponding eigenstate of the full HQET La- 
grangian (3). Although they correspond to the same energy eigenvalue (to the order 
in 1/m we are considering), they are not necessarily identical and in fact, it will be 
shown in the next section that they are related in the rest frame of the hadron by an 
unitary transformation. The normalization convention for the states \Mhqet{v)) is 
the usual one {MHQETip}\MHQETip')) = 2E^5^^\p ~ p'). 
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The matrix elements on the r.h.s. of (5) can be written with the help of the 
Gell-Mann-Low formula [0 as power series in 1/m. For example, the first term is 
equal to 

{Mhqet{v)\ I d?x{hHhH){x)\MHQET{v)) = {M{v)\ f d?x{hhi){x)\M{v)) 



+ -i{M{v)\T I d?x f d'yCooiy)Chihi){x)\M{v)) + ■■■ . (6) 



m 



Here, the field hi{x) is understood to be in the interaction picture where the term 
h{iv ■ D)h in the Lagrangian (3) is the unperturbed Lagrangian and the rest is an 
interaction. That is, we work to all orders in a^; this is motivated by our desire to 
have hadrons containing one heavy quark in the asymptotic initial and final states. 
Another useful splitting of the Lagrangian (3) is the one which keeps only the term 
hiv ■ dh as the free Lagrangian and treats the rest as a perturbation; this gives the 
usual double perturbative expansion in as and 1/m and is suited for studying the 
scattering of free quarks and gluons. Ci is the coefficient of 1/m in (3) and the 
states \M{v)) are eigenstates of the leading order term in the HQET Lagrangian. 
It can be seen that there are two distinct sources of mass dependence in (5)Q: 

a) the factor appearing explicitly in (5), which has a kinematic origin. It comes 
in through the very definition of a quantity of the theory, the number operator. 

b) powers of 1/m arising through application of the Gell-Mann-Low formula. These 
can be considered as having a dynamical origin, as they are due to the interaction 
terms in the Lagrangian (3). 

In the final expression obtained after application of the Gell-Mann-Low formula 
to (5), the mass dependence is concentrated into the factors of 1/m standing in front 
of matrix elements like those on the r.h.s. of (6). The latter are mass-independent 
quantities, properties of the HQET in the infinite mass limit. A similar expansion 
as a power series in 1/m with constant coefficients can be written for any other 
matrix element of heavy quark currents. In fact, it is this separation of the mass 
dependence one of the properties which makes the HQET useful. 

However, a closer look reveals some problems related to the form of the heavy 
quark number operator (4). Firstly, its expectation value in the unperturbed state 
\Miv)) 

{M{v)\ J d'xhj{x) \l + -^Pl] hj{x)\M{v)) (7) 



2mM 

differs from unity. Usually, the scale of a quantized field is set by requiring that the 
eigenvalues of the particle number operator be integer numbers. In our case, the 



^Neglecting the trivial m-dependence appearing through the normalization of the states. 
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matrix element (7) can be normalized to unity through the following field redefini- 
tion]]: 



h'{x) 



1 
2 



2m 



h{x) 



(8) 



Secondly, even more disturbing, the number operator (4) is not a constant of 
motion for the leading order HQET Lagrangian: 



dxo 



d xhj(x) 



1 + 



4m^ 



hjix)j^O. 



(9) 



To be sure, it is conserved when the dynamics is given by the full Lagrangian (3). 
The point is that we do not know how to solve the problem corresponding to this 
Lagrangian (nor can useful statements be made about it). We are forced to resort to 
the use of perturbation theory, which expresses any quantity of interest in terms of 
matrix elements in the theory defined by the first term in (3), the "free" theory. For 
these matrix elements, use can be made of the symmetries of this theory in order to 
extract useful predictions. 

Unfortunately, the nonconservation of the number operator (9) in the "free" 
theory has as a consequence the fact that the unperturbed states \M{v)) upon 
which the perturbative expansion in 1/m is built, simply do not exist. This means 
that the theory defined by (1-3) does not lend itself to a consistent approximation 
scheme in powers of 1/m. 

Fortunately, the solution to the first problem, the field redefinition (8), happens^ 
to be also the solution for the second one. In terms of the new field h'{x), the 
Eqs.(l-3) read 



Qix) 



+ 



4^2 



h'{x) 



h'l.h' + —h'{- P, 7^ + 7^ P,)h' 



1 ^ 



Am^ 



h' ( (t, pi 



Pi. V- Dh^- P, 7m 



h' 



iv- D 



-—Pi 
2m 



Am^ 



1 
2 



V- D 



'± +7m(^ pi -V- DP,] 
P, V DP, -\v- DP] 



(10) 

(11) 
h', 

h' . 
(12) 



This is precisely what one directly obtains by applying the method described in 
for deriving a heavy quark effective theory from QCD. The number operator (5) is 

similar transformation is sometimes made when deriving the PauU equation to order \/m?. 
See e.g. |. 

■^That is, the same is not true at higher orders in 1/m. 
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written now simply as 



7V = y d^x{h'h'){x) 



(13) 



and is a constant of motion for the free theory (i.e., when keeping only the leading 
term in (12)). 

Equally important, it stays conserved up to each order in the perturbative expan- 
sion in 1/m given by the Lagrangian (12). As discussed above, this is a prerequisite 
for the consistency of the approximation scheme. To see this to order 1/m, we note 
that the equation of motion for the field h'{x) can be written as 



IV' 



1 
2m 



h\x) = . 



This can be used to write 



-N 



d^x(h' D -vh' + h'v Dh') (x) 



(14) 



(15) 



2m 



I d^x (h' pl h' - h' pl h'j (x) = , 



where an integration by parts has been performed, followed by the neglect of the 
surface terms at infinity. 

The number operator for heavy quarks at rest is conserved also at order 1/m^ in 
the effective Lagrangian (12). This can be seen most easily by rewriting the 1/m^ 
term in the latter as 



^ P1 [v D, Pj] + V- D] I ( P^ v^YF,. - v>^rF,. P, 



This gives an equation of motion for h'{x) 



^ 2m 8m^ 



P± v^YF^^u - v^YF^u P_ 



h\x) = 0, 



(16) 



;i7) 



which does not contain, in the rest frame of v, any time-derivative acting on h'{x), 
except in the leading-order term. The above procedure can be repeated therefore, 
step by step, for this case as well. 

Everything which was said up to now refers to the number operator for heavy 
quarks at rest. However, it can be seen that similar results hold for the number op- 
erator for heavy quarks moving with an arbitrary velocity v. Heavy quark effective 
theory (even truncated to a finite order in 1/m) is invariant under Lorentz transfor- 
mations, provided the velocity v is transformed at the same time. Furthermore the 
number operator is a Lorentz invariant and is the same in any reference frame. In 
particular, if it is conserved in one reference frame, it should also be conserved in 
any other reference frame. 
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3. Comparing matrix elements of the currents 

We have shown in the preceding Section that the HQET (1-3) cannot be formulated 
as a consistent perturbative expansion in powers of 1/m. In principle, this would 
suffice to reject it. However, we would like to present in this Section a different, more 
pragmatic point of view, which will permit a better understanding of the difference 
between the two HQETs, making at the same time the connection with the, perhaps 
more familiar, theory of the point field transformations in the Lagrangian formalism. 

We will compare in the following the expHcit predictions of the two HQETs, 
given by (1-3) and respectively by (10-12), for the matrix element of the current 
cTb. Here b, c are two different heavy quarks moving with velocities v, v' and having 
masses mbjiric. By examining the form of the currents and Lagrangians, it is clear 
that a possible difference can only show up at order (9(l/m^). This is equal to 

{M'HQET{v')\{crbm\MHQET{v)) iW-12) - (M^q^^(^;') | (cEfe) (0) |M^,qet(^;)) (1-3) = 

^,{M'{v')\T J d'xh^:}{x)iP'l v'. D +v'- B^i)/i5^(x)(;iS^r/i(''))(0)|M(^;)) 

+ ^(M'(^')|T / d'xh^:\x){Pl V D +v- DPl)h^\x){h^)Th^^mmv)) (18) 

-^(M'(t;')l(/^S^ Pi Th^^m\M{v)) - ^^{M\v')\{h^)T Pi h^^'m\M{v)) . 
We can rewrite this as 

"8^^^'^^'^'^ / (^^'^ ^ + ^ ^^'^) r^i'^)(o)|M(t.)) 
-^(M'(^')l(^? k rh^:m\M{v)) - ^mv'm^^r pl /.l^))(o)|M(^)) , 

(19) 

where C denotes the sum of the "free" Lagrangians of the two heavy quarks: 

C = h^^}iv' ■ Dh^) + U^Hv ■ Dh^'^ . (20) 

The matrix elements (19) between hadronic states can be expressed with the help of 
the reduction formalism as the residues at the corresponding bound-state poles of 
the appropriate Green functions. For example, the first term in (19) can be obtained 
from the Green function 

(oiT(g-/.j))(y) / d'x [h^) p'i^) + ^)Pi {x){h^:)vh^))mm){m 

(21) 
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where [qh^^)) and (h^^^'q) are the interpolating fields of the outgoing and respectively, 
incoming hadrons. Only their fiavor structure has been made explicit, the Lorentz 
and Dirac structure needed to endow them with the correct quantum numbers of 
the respective bound states are implicitly assumed. One can apply now the relation^ 

(0|T (x)gi(xi)g2(x2) . . . |0) = (22) 

t Y: 5^'\x - x.){0\^Q,ix^) . . . I^F^^ (x,) . . . |0) , 

to transform (21) into 

t{0\T{qh^}){y){h'^:> P'l Th^^)mh^^^q){m (23) 

+z(0|T(g P'l h^:^){y){h^)vh^^)mh^\){m ■ 

The first term can easily be translated back into a matrix element between hadronic 
states, equal to 

z{M\v')\{h^:^ p'lvh^^mmv)) . (24) 

In a completely analogous way, the second term in (19) can be transformed to the 
form 

^(0|T(g/.5))(y)(^l?)r^W)(0)(^(^) Pi q){zm (25) 
+z(0|T(g-/.S))(y)(/.J)r Pi hf)mh^h){zm , 

of which the second term gives again a matrix element contributing in (18) a quantity 

i{M\v')\{h^:)v Plh^\Q)\M{v)) . (26) 

Collecting now everything, one can see that the terms (24) and (26) exactly cancel 
the matrix elements of the two local terms in (18). The difference (18) of the matrix 
elements of the current cTh in the two HQETs is proportional to the residue at the 
corresponding bound-state poles of the Green function 



^(0|T(g P'l h^}){yWTh'^^)mh^^h){zm (27) 
8m^ 



Given in |^ and [ pO[ with a wrong coefficient on the right-hand side. The correct form can be 
found for example in . 
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Next we observe that 

Plh^:> = {d' - W ■ Df + ^-a.^F^'^ + igv'p^F^'^^ /.^^ . (28) 

If we keep in this relation the terms containing only the heavy quark fieldQ, the new 
interpolating field {qh^^)) {qd\h^^)) has the same quantum numbers as the old 
one and is therefore equally qualified to represent the same hadronic bound states 
as the latter. A similar argument can be applied to the second term in (27). We are 
led therefore to conclude that the two Green functions in (27) have nonvanishing 
residues at the bound-state poles associated with the hadronic states \M{v)) and 
\M'{v')). Hence the two HQETs give different answers to order 0{l/m^). 

In taking the difference (18), the matrix elements on the r.h.s. are taken between 
the same states. These are the physical states in the infinite-mass limit. On physical 
grounds these states exist and are unique. The reason for obtaining different results 
with the two HQETs should be clear: in passing from one theory to another, we 
have "forgotten" to make the change of variable (8) in the interpolating field for 
the mesons too. Once this is done, the contribution from the second term in (8) 
will cancel the matrix elements (27). But this means that the two interpolating 
fields {qh) and {qh') cannot be simultaneously good interpolating fields for the same 
states, to be used with their respective theories. That is, we have an alternative: i) 
Use the HQET (1-3) with the interpolating field [qh) for the heavy meson. Then, 
when working with the HQET (10-12), the correct interpolating field to be used 
should be 

m-^,mih'). (29) 

8m^ 

Or, ii) Use the HQET (10-12) with the interpolating field {qh') and the HQET (1-3) 
with the interpolating field 

m + ^.mih). (30) 

8m^ 

Only one of these possibilities can be true (or none). We will prove in the next 
section that the alternative ii) is the correct one, by showing that the normalization 
of the interpolating field {qh') is the same as that of its QCD counterpart {qQ): 

{mh'H)\MHQET) = {0\{qQ)\MQCD) . (31) 

The situation can be understood of inapplicability of the well-known 

equivalence theorem According to this theorem, any field transformation 0' = 
(j) + F{(j)), where F{(j)) contains only terms with at least two fields, leaves the physical 

^The terms containing at least two fields do not give a pole and vanish under multiplication 
with the LSZ factors. 
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content of the theory (the S-matrix) unchanged. In other words, any arbitrary field 
(f)' is as good as (j) as an interpolating field, as long as it has the same normalization: 

(O|(/.'(a;)|0(p)) = (O|0(x)|0(p)). (32) 

The condition on the field transformation -F(0) has precisely the function of assuring 
that the new field is correctly normalized. 

In our case, the two theories are related by the field transformation (8). It is 
easy to see that the above-mentioned condition on the field transformation is not 
satisfied for this case and therefore the two fields h{x) and h'{x) do not have the 
same normalization. We will show now that it is the field h'{x) which is correctly 
normalized. 

4. Proof to all orders in 1/m 

In this Section, after a short review of the results obtained in 0, we will take up 
the arguments of the preceding two Sections and will prove that they are true to all 
orders in 1/m. 

In the approach of 0, the heavy quark field h{x) in the effective theory is 
expressed in terms of the QCD field Q{x) as 

h{x) = h^^\x) + h^-\x) (33) 

with 

/.W(x) = ^-^e'^--^{...e-^^^e-^^^e-^^^)Q{x) (34) 

h^-\x) = l^e-^"^^-^(---e-2^^3''e-2^^2^e-2^^i')g(x). (35) 

The differential operators Of are determined so that the Dirac Lagrangian 

C-Dirac = Q{ip - m)Q , (36) 

does not contain, when expressed in terms of the fields h'^^\x), any terms which 
couple /i^"*"-* {x) with h'^~'> {x) . This is done by removing order by order in 1/m, through 
succesive field transformations, every term in the Lagrangian which anticommutes 
with i/. These represent the generalization of the "odd" operators in the usual 



Foldy-Wouthuysen procedure. The first few Of^ are 

Of = iP^ (37) 

Of = -\{P^P\\ + P\\P^) (38) 

Of = -i{\pl + \p^P\ + \p\\P^P\\ + \p\P^). (39) 
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All of satisfy {Of^ ^ = and contain, in the rest frame of only spatial deriva- 
tives which act on the Q field in (34,35). 

It can be seen that this is a generalization of the usual Foldy-Wouthuysen trans- 
formation which decouples the upper two components of the Dirac spinor from the 
lower two ones. In a similar way, the transformation (33-35) decouples those com- 
ponents of the Dirac spinor which correspond to a particle moving with velocity v 
from its components corresponding to an antiparticle moving with the same velocity. 
In the rest frame of v the usual Foldy-Wouthuysen transformation is recovered and 
the equation of motion for the Heisenberg field h^'^^[x) becomes simply the Pauli 
equation. Motivated by the usual physical interpretation of the Foldy-Wouthuysen 
procedure, one can say in a loose way that heavy quark effective theory represents 
a "non-relativistic approximation about some arbitrary velocity f " . This has been 
formulated in a precise language to order l/m*^ by Grinstein [1^ (in the rest frame 
of v). 

The recursive algorithm for determining the operators Of and thereby the La- 
grangian of the effective theory given in is very simple and lends itself easily to 
implementation in a symbolic manipulation language (e.g. FORM [|14|)- The effec- 
tive Lagrangian obtained in this way to order (9(l/m^) reads (only for the particle 
field h^^^ix), which will be written simply as /i(x)) 

^ = E -7 A (40) 



i=0 



with 



£0 = hiv-Dh (41) 



A = -h(^D'-{vDy + ^-a,,F^'')h (42) 



3 



;h{D^-{v-Df + ^-a^^''F^,^ h. (44) 



The equation of motion for the h{x) field in the rest frame of the heavy quark is 



= {gA'H^ + 7^"'"^' - ^cr . B^t^ - -^(divE» + gU^A' 
ot 2m 2m 8m^ 



• 2 

■ lotm'' - ^fabc^ ■ (A" X - /-a ■ X (P - gA't') 

im^ Qm^ 4m^ 



8m^ 



(P - gA^ey - ga ■ B^f 



+ 



8m^ 



)h 



+0(l/m*) , (45) 



10 



with P = -zV, F;, = d^A: - d,A; - gfabcA^Al and D^F^^^ = d^F-, - gfabcA\F^^,. 
The electric and magnetic color fields are defined as E''"- = 5*"^ = —^e^^'^Fj-j^. 
This agrees with the known form of the Pauli equation 0, (though, to our 
knowledge, the complete 0{l/m^) contribution is new) although the effort involved 
in deriving it, even in QCD, is appreciably smaller than the one required by the 
traditional method using unitary transformations. In the Appendix are given the 
operators which allow the determination of the HQET Lagrangian and of the 
Pauli equation up to 0(l/m^). The corresponding expressions have been calculated 



up to 0{l/m ) and will be presented elsewhere ^7 



In Sect. 2 we have shown that the number operator in the HQET of has 
the simple form shown in Eq. (13). An important consequence of this fact is that 
its expectation value in the unperturbed state \M{v)) is equal to 1. Actually, this 
property is preserved to all orders in 1/m. To see this, we note from (33-35) that 
the number operator can be written as 



N = J d-'xQ^oQ = (46) 
'^i^h^imv ■ x^^^-tmv -xi/j^^ f (/,(+)/,(+) _ /,(-)/,(-)) (a;) . 



Here we have made use of the above-mentioned property of the differential operators 
of of containing only spatial derivatives in order to integrate by parts and drop the 
surface term. In the second step the relation {Of', = has been repeatedly used 
(remember i/ = 70). This proves the desired result. 

Another crucial property of the HQET in which was proven at order l/m^ 
in Sect. 2 is the fact that the number operator in the full theory is a constant of 
motion order by order in 1/m. By making use of the preceding result concerning 
the form of the number operator, one can see that this is indeed the case to any 
order in 1/m. The proof for this is similar to the one given in Eq. (15), where the 
main ingredient was the absence of the time derivatives in the HQET Lagrangian 
(in the rest frame of v). It has been proved in 0] that the method of constructing 
the effective Lagrangian does not introduce generally any v ■ d derivatives on the 
heavy quark field. Any extra number of derivatives can be reexpressed in terms of 
F^jy and its derivatives. Thus, the conservation law of is guaranteed to hold at 
any order in 1/m. 

The conservation of the number operator order by order in 1/m has one interest- 
ing consequence, in that its expectation value remains equal to unity order by order 
in 1/m. This follows from the nonrenormalization theorem for conserved quantities, 
according to which matrix elements of conserved operators have the same values as 
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in the noninteracting theory. We have just proved that the expectation value in the 
"free" theory is unity, so that this is also the value to each order in 1/m. This is, 
of course, what one expects from the complete theory (QCD), and is reassuring for 
the consistency of our scheme. 

The sequence of transformations (33-35) can be regarded as the effect of an 
unitary transformation on the Hilbert space of the initial (QCD) theory. Again, we 
consider the rest frame of v. The respective unitary transformation is, of course, the 
usual Foldy-Wouthuysen transformation in second-quantized form. We have 



\Mhqet) = U\Mqcd) , h = /iW + /i(-) = UQU"^ , (47) 



with 



U = exp \-imxQ{N^+^ - Ar(-))1 . . - exp [-^ / d?xQ'^O^Q'\ exp [-^ / d^xQ'^O^Q 

\_Ziflh J J \_Ziflh J 



Here 



i2m J O y]e 2m J ^ e 2m 



Q'{y) = e2m J ^ Q{y)e 2m J ^ ' ^ ^ e 2m ' Q{y) , (49) 

Q {y)=e2m^J Q {y)e 2m^ J =e 2m^ Q'{y),{50) 

etc. , 

ft<+>® + A<->(?)= (51) 

give the effect of the successive transformations on the heavy quark field operator 
and 

7V(±) = I d^x/iW/iW (52) 

are the number operators for heavy particles and antiparticles, respectively. In 
these relations all the operators are taken at equal times and = Qjo- The second 
equahty in Eqs. (49-51) follows through the use of the canonical anticommutation 
relations for the heavy quark field. Here an important point is that these relations 
arc left invariant by the succesive unitary transformations. 

We can see now that the relation anticipated in the previous Section 

{0\{qQ)\MQCD) = mqh)\MHQET) (53) 

holds indeed true. This shows that the field h{x) defined by the sequence of trans- 
formations (49-51) has the same normahzation as the QCD field Q{x) and therefore 



.(48) 
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the results obtained with the HQET formulated in terms of it are the correct ones. 
On the other hand, as already remarked, the field corresponding to the usual HQET 
(1-3) has a different normalization. 

One curious property of the effective Lagrangian (40) is the fact that the number 
operators A^^"*") and A^*^") for heavy quarks and respectively, antiquarks, are sepa- 
rately conserved to any order in 1/m. This can be traced back to the invariance of 
the Lagrangian under the field transformations h —>■ exp{ia)h and h —>■ exp{i]/a)h, 
the generators of which are related to the two number operators. But this is surpris- 
ing since we know that in the full theory (QCD), it is only the difference of these two 
operators which is conserved (electric charge conservation), but not each of them. 
That is, HQET has one more integral of motion than the underlying theory. This 
seems to indicate the fact that the predictions of HQET do not approach arbitrar- 
ily close those of QCD when infinitely many orders in 1/m are taken into account. 
Namely, effects connected with heavy quark pair production are not described, to 
all orders in 1/m. See also Ref. [|18|. We mention that this limitation applies only to 



the tree-level HQET Lagrangian and it is possible that matching corrections might 
remedy this deficiency. If not, this would signal a limitation in the predictive power 
of HQET for very large orders in 1/m. 

5.Renormalization of the order Lagrangian 

There is a class of contributions at order (9(l/m^) which are not described by the 
HQET Lagrangian (43). These arise in the full theory (QCD) from graphs contain- 
ing one heavy quark in closed loops. They are therefore quantum effects which have 
not been accounted for in the tree level (classical) effective Lagrangian (40). Con- 
sequently, the missing contributions must be added by hand as gluonic operators 
with coefficients which have to be determined from a comparison (matching) with 
the full QCD graphs. For simplicity we treat only the case when no light quarks are 
present. There are three possible operators which can appear when inserted at zero 
total momentum: 



Oo = -{D^F;,){D^F^''^)--g{D^F;,){hYm (54) 

O2 = \gfa,cF;^,Fl^F;^. (56) 

The operator Oq vanishes when the equations of motion for the gluon field are used. 
By making use of the Bianchi identity 

D^F;^ + D,F:, + D^Fl = (57) 
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one can see that the following relation between these operators holds true, when 
inserted at zero momentum: 



0, + 202 + l{D>^Fp{D,F 



a\v\ 



0. 



(58) 



The HQET Lagrangian (43) supplemented by these terms reads 



m 



(59) 



Here we have denoted by ^3^4 the spin-symmetry violating terms in (43): 



O3 
O4 



(60) 
(61) 



Additional 0(l/m^) contributions also appear from double insertions of 0{l/m) 
terms in the effective Lagrangian: 



O5 



J d^xTi^ha ■ Fh){x){^ha ■ Fh){0) 



- I d''xT{h{iDfh){x){h{tDfh){0) 



The coefficients Ci(/i) obey a renormalization group equation 



(62) 
(63) 
(64) 



/^^Ci(/i) +7jiCj(/i) = 0. 
The anomalous dimension matrix with z, j = 1 — 7 has the form 



7 



167r2 



/ -22 -16 \ 

9 4 

-12 

12 

48 48 -6 

-12 

y oooooy 



(65) 



(66) 



The nonlocal operators ^5,6, 7 are not required by the local ones as counterterms 
and therefore their coefficients 05^6,7 (/^) evolve independently from the latter. Their 
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anomalous dimensions (the lower 3x3 diagonal block) have been calculated previ- 
ously IC]. The renormalization of the dimension-6 gluon operator O2 has been 



studied in the presence of a light quark in connection with the renormal 



ization of the corresponding condensate appearing in QCD sum rules. We agree 
with the results of these works for the case of a massless quark. However, when 
considering a HQET heavy quark, the renormalization is different. The diagrams 
to be calculated for the renormalization of Oi^2 are shown in Fig.l, those for 03^4 
in Fig.2 and in Fig. 3. The necessary algebraic manipulations have been done with 



the help of FORM |1J] . We have used for the gluon field a Feynman version of the 



background field gauge [^, where the combination gA^ is not renormalized. In 
general, some other operators will be needed as counterterms, which vanish when 
the equation of motion for the heavy quark is used: 

Os = '-lha-F{vD)h~'-^h{D-v)a-Fh (67) 

Og = ihD^{v D)h - ih{D ■v)D'^h (68) 
Cio = h{iv-Dfh. (69) 

Their contributions have to be properly separated. Even if these operators do not 
contribute to the hyperfine splittings, they do appear when considering corrections 
to the matrix elements of currents. Therefore, a complete study of the 0{as/'m?') 
order corrections should take them into account. 

The problem of the renormalization of the 0{l/m?') effective Lagrangian has 
been treated previously by Lee p. Our results do not agree with his. For example, 
we find that the operators 03,4 are renormalized always in the combination O^ + Oi^, 
whereas in they appear together as 2(^3 + (denoted as 02)- Unfortunately, 
this combination is not even hermitian (whereas + is) . Its expectation values 
on the other hand have a direct physical interpretation in terms of hyperfine level 
splittings which are clearly real quantities. 

The coefficients q at the matching scale n = m are 

Co(m) = --, ci(m) = -- + 



4' 4 ' 1207r2 

C2(m) = + C3(m) = l (70) 

C4(m) = 1, C5(m) = C6(m) = ^. 

Here the leading terms of the coefficients co,i.2(^) have been determined by re- 
expressing the first term in (43) in terms of Oo.i,2 by making use of (58). The 
contributions proportional to in Ci^2{^) constitute the quantum corrections re- 
ferred to at the beginning of this section. They have been obtained by computing 
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the QCD graphs in Fig. 4, where the fermion in the loop is a dynamical one of mass 
m. 

To one-loop order, the combinations of Oi^2 which are multiplicatively renormal- 
ized are 

C»'i = O1 + 2O2 (71) 
O'^ = 9C1 + 8C2 (72) 

with the anomalous dimensions 7^ = —^,72 = "If? and the initial conditions 

c[{m) = — I + 2^QQ^2 ,c'2{m) = —7200^' "^^^ matrix elements of O'l can be related 
via (58) to those of the first term in (43). Therefore, as far as its matrix elements 
are concerned, the latter operator can be considered as being multiplicatively renor- 
malizable with the anomalous dimension 7^ Summing up the leading logarithms in 
c[ 2(/^) and C3,4(/i) with the help of the renormalization group equation (65) and the 
initial conditions (70) gives the effective Lagrangian £2 for an arbitrary n < m: 

£2 = -^(^y^o'2+(^y"'-{hv^t^h)D^F;:^ 



+ 



ISOOvr \as{m) 
' ajjjy 



_3_ 

' 11 




(73) 



Here we have used the running of the coupling proper to a pure gauge field. 

It is interesting to observe that the first term is related to what we would have 
obtained if we had "completely" integrated out the heavy quark. In the absence of 
any other (light) quarks, there is only one independent dimension-6 operator which 
can appear into an Euler-Heisenberg Lagrangian for the gluon field, which can be 
chosen to be O2. Since now the gluon field is a free one, Eq. (58) gives O2 = — IOC2 
0. Substituting this into (73) yields the well-known result ||2^ 



Also, 72 is identical with the anomalous dimension of O2 for a pure gauge field ||22 

However, keeping the first term in (73) is an inconsistent procedure, as there are 
contributions of comparable magnitude into the second term which were neglected. 
They are the non-leading logarithms in c[ of the form a<j(Q;s ln(/i/m))", whose sum- 
mation by the renormalization group will require the knowledge of the anomalous 
dimension matrix (66) to two-loop order. This is somewhat a disappointing feature, 
since, in the view of the preceding paragraph, they constitute the genuine vacuum 

^This equality is to be understood in the sense that the two operators have identical matrix 
elements. 
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polarization contribution to the heavy quark effective Lagrangian. Nevertheless, 
such a two-loop calculation is beyond the scope of this paper. 

Unfortunately, there are reasons to expect that the results for the anomalous 
dimensions (66) will be changed when light quarks are included. There are 19 new 
dimension-6 operators including both light- and heavy-quark fields and even purely 
light-quark operators which can be induced by renormalization in the presence of 3 
massless quark species. Thus, there are 6 operators of the form 



with the two light quark fields qq, q = {u,d,s) in a SU(3)-singlet combination. 
There are 12 4-light-quark operators of the following two types 



with two possibilities for each of them, corresponding to the two occurences of 1 in 
the product of representations (1©8)®(1©8). Here T = (1 , 7^ , cr^j,). And finally, 
1 operator of the type iq{F^^ ,'-j^D^}q. For a similar calculation in QED with one 
light and one heavy fermion, see Therefore the results of this section are mainly 
of academic interest. 

6. Applications 

In this section we discuss some of the implications of our results. One possible 
application would be to the hyperfine splittings of hadrons containing one heavy 
quark. As is well-known, in the infinite quark-mass limit mq 00, the states of 
such hadrons can be labelled by simultaneously specifying the angular momentum 
and parity of its light degrees of freedom and the total angular momentum 
S = sq + St (with sq the heavy-quark spin). Since in this limit the spin of the 
heavy quark decouples, there exists a degeneracy between the two states with the 
same but different S = si± \ (for st 7^ 0). The effect of a finite heavy quark 
mass is to lift this degeneracy, leaving us with the observed hadronic spectrum 
consisting of closely spaced doublets. The splittings within these doublets are mainly 
a l/mg-effect due to the spin-symmetry violating term in the effective Lagrangian 
Ci (42). At order l/wig, an additional contribution is expected to appear from the 
corresponding terms into C2 containing a^j^p. We can get some insight into the size of 
these contributions by examining the role played by these terms in the constituent 
quark model. In the rest frame of the heavy-quark, they can be seen to be responsible 
for the well-known spin-orbit interaction between the spin of the heavy quark and 



{qq){hh) , {q'j^q){hv^h) , {qa ^„q){ha ^^h) 
{qt''q)Cht''h) , {q-f^t''q){hv^t''h), {qa^Xq){ha^Xh) 



(75) 



{qVq){qVq), {qTt^q){qTt''q) 



(76) 
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its orbital angular momentum. Writing the spin- dependent part of the Breit-Fermi 
Hamiltonian for a pair of a heavy- and light-quarks as ^6 



with 

^spin-spin ^ ^!^5(3)^^)^ . ? 

^tw^A '' '~ — w — J ^ ^ 

1 1 

-r X ■ Sg -t- — f xpg-SQ, (79) 



then it is the first term in ffspm-orbit -y^^j^ic]^ arises from the above-mentioned terms in 
£2- Here eQ^rriQ^eq^mq are respectively, the heavy and light quark electric charges 
and masses and is — | in mesons and — § in baryons. The constituent quark model 
suggests thus that the £2^contribution to the mass splittings of s-wave mesons hke 
rriD* —rriD or ruB* —rriB is expected to be very small, being exactly zero in the valence 
quark approximation!]. This provides an explanation for the fact that the observed 
hyperfine mass splittings of these mesons are well described by a simple l/mg law 
(including the leading order logarithmic corrections) 0. It should be manifest on 
the other hand in the corresponding hyperfine splittings of the p-wave mesons which 
therefore might display larger deviations from a simple l/mq law. Even for this case, 
the spin-spin forces will probably give the dominant contribution. 

The £2-terms can be expected to play an important role for the case of the first 
excited negative-parity 1 = baryons with sj*^ = 1~. In a constituent quark model 
P5| , these states are described as having the light quark pair (m, d) into a spatially 
symetrical state and their total spin is zero. If the inter-quark forces are taken to be 
of a harmonic oscillator type, the relative angular momentum of the two light quarks 
is zero and the resulting diquark "looks" like a scalar particle which orbits about the 
heavy quark in a p-wave. When the heavy quark spin is added, two baryon states 

1 3 

will result Aa" and Aa". The only term at 0{l/mQ) responsible for the hyperfine 
splitting between these two states is the last one in ^«p«"-o'"''«* (79)|. One may thus 

^Of course, there are also 1/m^-order contributions arising in second-order perturbation theory. 
If one takes the usual point of view according to which the Breit-Fermi Hamiltonian is only to be 
used as a first-order perturbation, they vanish (for a review see |p7|| ). 

®This statement is not specific to the harmonic force case. It is only based on the symmetry of 
the wave-function under a permutation of the two light quarks. 
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expect the £2^erm (the first one in (79)) to contribute an appreciable fraction to 
the hyperfine splitting of these baryons. 

Another application of our formalism is connected to the possibility of a model- 
independent extraction of the Kobayashi-Maskawa matrix element Vet as described 
in [^. This method is based on the absence of l/rric corrections to the form- 
factor /iAi(l) describing the matrix element of the axial current = 07^756 for the 
transition B —>■ D*eu: 

{D*{v\e')\c^,^,b\B{v))= (80) 
yJuiBmn* [/lAi {v-v'){l + v- v')e'* - Ha^ {v ■ v') {e'* ■ v)v^ - Ha^ {v ■ v') (e'* ■ v)v'^ . 

In the infinite-mass limit, hA^{l) = 1 and the corrections to this result appear only 
at order 1/m^ It has been argued in that these corrections amount to about 
2%. As mentioned in the Introduction, the calculation of is based upon the 
theory defined by (1) and (3). We have seen that these relations are incorrect and 
have to be replaced by (10) and, respectively (12). These changes imply a different 
form for the axial current in the HQET which contains now two new contributions: 

A^. = -^h^) l,l.h^^ - 3^^S^.75 P^' , (81) 

proportional to 1/m^ and respectively 1/m^. It is natural to ask how this change 
will affect the predictions of It turns out that, due to the particular way of 
evaluating the 1/m? corrections used in they remain unchanged. A simple way 
to see this is presented in the following. For notations we refer to 0]. 

The form-factor /iyii(l) at the equal-velocity point f ■ f ' = 1 is equal to 

hA,{l) = 1 + €%{!) + eliiil) + e,efe(m2(l) + mg{l)) + 0{e') , (82) 

with eg = l/(2mQ) and £1,2, rn2,9 are functions of v ■ v'. Vector current conservation 
gives two conditions on the values of these functions at the equal-velocity point 
V ■ v' = 1: 

2£i(l) + (mi(l)-m8(l)) = (83) 
2^2(1) + (m4(l) - mn(l)) + {m,{l) - mi2(l)) = . (84) 

The (9(l/m^) correction to /iAi(l) is 

/lAi(l) - 1 = (ec - e6)M2(l) - efe^i(l)] + e^e^A (85) 

with 

A = ^i(l)+^2(l)+m2(l)+m9(l). (86) 
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The constants ^1,2(1) are determined from a comparison with the predictions of the 
Isgur, Scora, Grinstein and Wise model. This determination is independent of the 
precise form of the HQET used. A can be rewritten by making use of the relations 
(83-84) as 

A = -^(mi(l)-m8(l))-^(m4(l)-mn(l))-^(m5(l)-mi2(l)) + (m2(l)+m9(l)) . 

(87) 

Now it can be observed that the correction (81) only affects the internal structure 
of the £i_6 form-factors (see (4.34) in ||^) and has no effect on the structure of the 
functions mi_24. Therefore the prediction for the value of the constant A resulting 
from (87) remains valid. 

However, any calculation based on first principles of i'l-e? such as with the help 
of QCD sum rules or on lattice, should give an incorrect result when the HQET 
(1-3) is being used. 
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Appendix 

We give here the differential operators O^^, needed for constructing the HQET 
Lagrangian up to 0{l/m^). They are 



5 ^ 3 ^ ^ ^9 3 ^9 ^ ^ 1 

16^ 



Ot = ^{\P\ + ^PIP\ + \PIP\\P^P\\ 

+ ^PIPIP^ + ^P^P\\PIP\\ + \P^P\\P^P\\P^ + ^p^p\pl 



+ I^P^PI + \P\\PIP\\ + ^P\\PIP\\P^ + \P\\P^P\\PI 
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Figure Captions 



Fig.l Diagrams contributing to the anomalous dimension matrix of Oi and 02- (a) 
and (b) have a symmetry factor 1/2 and (f) appears together with another diagram 
obtained by reversing the direction of the fermion hne. In (d-g) the blob can only 
be an insertion of Oi. 

Fig. 2 Inserting O3 or O4 in heavy quark vertex functions. Diagram (e) has a 
symmetry factor 1/2. Three other diagrams contribute which are obtained from 
(a,b,d) by reversing the direction of the fermion line. 

Fig.3 Double insertions of Ci. The square and the circle may represent either 
different or identical terms in (42). To each diagram which is asymmetrical with 
respect to its central vertical axis, its mirror image should be added. 

Fig.4 Diagrams needed for the matching conditions (70). 
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